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We propose a new protocol for the manipulation of a three-level artificial atom in Lambda (A) 
configuration. It allows faithful, selective and robust population transfer analogous to stimulated 
Raman adiabatic passage (A-STIRAP), in last-generation superconducting artificial atoms, where 
protection from noise implies the absence of a direct pump coupling. It combines the use of a 
two-photon pump pulse with suitable advanced control, operated by a slow modulation of the phase 
of the external fields, leveraging on the stability of semiclassical microwave drives. This protocol 
is a building block for manipulation of microwave photons in complex quantum architectures. Its 
demonstration would be a benchmark for the implementation of a class of multilevel advanced 
control procedures for quantum computation and microwave quantum photonics in systems based 
on artihcial atoms. 

PACS numbers: 03.67.Lx,85.25.-j,42.50.Gy 


Advanced control of multilevel quantum systems is 
a key requirement of quantum technologies [1], en¬ 
abling tasks like multiqubit or multistate device pro¬ 
cessing [2-4] by adiabatic protocols, topologically pro¬ 
tected computation [5] or communication in distributed 
quantum networks [6-8]. These are currently investigated 
roadmaps towards the design of fault tolerant hardware, 
i.e. complex quantum architectures minimizing effects 
of decoherence [9, 10]. In this scenario artificial atoms 
(AAs) are very promising since, compared to their natu¬ 
ral counterparts, they allow for a larger degree of integra¬ 
tion [11-14], on-chip tunability, stronger couplings [15] 
and easier production and detection of signals in the 
novel regime of microwave quantum photonics [16]. De¬ 
coherence due to strong coupling to the solid-state en¬ 
vironment [6] is their major drawback. Over the years 
it has considerably softened [10] yielding last-generation 
superconducting devices with decoherence times in the 
range ~ 1 — 100 ^s [2, 4, 20]. 

Gombining potential advantages of AAs is by no means 
straightforward. Protection from decoherence often im¬ 
plies strong constraints to available external control, 
which pose key challenges when larger architectures are 
considered [14]. In this work we study a simple and 
paradigmatic example, namely a three-level AA driven by 
a two-tone electric field in the Lambda (A) configuration 
[Fig.1(a)]. Implementation of this control scheme in last- 
generation superconducting hardware may in principle 
benefit from low decoherence, which however is achieved 
at the expenses of suppressing the direct coupling of the 
pump field, and of possible limitations of selectivity in ad¬ 
dressing specific transitions. In this work we show how 
to implement an efficient A configuration in these con¬ 


ditions, and we propose a dynamical scheme allowing to 
operate quantum control. This solves the problem raised 
in the last decade by several theoretical proposals on the 
implementation of advanced control by a A-scheme in 
AAs [3, 21-25], which still awaits experimental demon¬ 
stration. 

Quantum control via a dynamical A scheme is very 
important because it may provide a fundamental build¬ 
ing block for processing in complex architectures. In¬ 
deed adiabatic evolution may be used to trigger two- 
photon absorption-emission pumping cycles, which al¬ 
low for on demand manipulation of individual photons 
in distributed quantum networks, as proposed in the 
cavity-QED realm [26, 27]. Demonstrating control by 
a A configuration in last-generation AAs would extend 
this scenario to the microwave arena, opening the per¬ 
spective of performing demanding protocols in highly in¬ 
tegrated solid-state quantum architectures [13, 14], which 
are usually subject to specific design constraints [28]. 
Examples are adiabatic holonomic quantum computa¬ 
tion [29], information transfer and entanglement genera¬ 
tion [24, 30, 31] between remote nodes, and other sophis¬ 
ticated control protocols [32]. 

The A scheme is described by the standard Hamilto¬ 
nian in the rotating-wave approximation (RWA), which 
in a double rotating frame reads [9] 

H = ^|0)(21 + ^11)(21 +h.c. + <jll)(ll +<5,12)(21 . (1) 

Here Hp g are the Rabi frequencies of the (pump and 
Stokes, respectively) external fields, with single-photon 
detunings Sp^s, and S = 6p — 5s is the two-photon de¬ 
tuning [Fig. 1(a)]. This latter is a key parameter, since 
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for (5 = 0 the system admits an exact dark state \D) := 
+ (ris|0) — rip|l)): the system is trapped in \D) 

despite the two fields triggering transitions towards |2), 
a striking destructive quantum interference phenomenon 
named coherent population trapping [9, 34]. Sensitiv¬ 
ity to 5 is critical since no exact dark state exists for 
(5^0, even if partial trapping is still supported. Con¬ 
trolling the dynamics of \D) would lead to the observa¬ 
tion of basic interference effects, allowing also important 
applications. A representative example is stimulated Ra¬ 
man adiabatic passage (STIRAP) [9, 35], a powerful tech¬ 
nique allowing faithful and selective population transfer 
in atomic physics. By adiabatically varying in the 

so called counterintuitive sequence [Fig.1(b)], the state 
\D{t)) evolves from jO) to Jl), in the absence of a direct 
coupling and never populating the intermediate state j2). 
STIRAP is a benchmark for multilevel advanced control. 
Its robustness against imperfections and disorder may al¬ 
low to develop new protocols [2, 32] with important ap¬ 
plications in hybrid networks, composed of many AAs or 
microscopic spins interacting with quantized modes [36]. 

Several works in the last decade [3, 21-25] proposed 
implementations of A-STIRAP in superconducting AAs. 
However dynamics in the A scheme has not yet been 
experimentally demonstrated for a fundamental reason: 
protection against low-frequency noise [6] necessary to 
achieve large decoherence times is attained by operating 
the device by a Hamiltonian with exact or approximate 
symmetries [5, 37]. In particular low decoherence in last- 
generation superconducting AAs is obtained by enforcing 
parity symmetry, which however implies cancellation of 
the coupling to the pump field [3, 21, 23, 25, 39]. A- 
STIRAP could be observed by breaking the symmetry of 
the device [21, 23], but at the expenses of an increased 
noise level. Analysis of a case study [7] has shown that 
efficiency, i.e. the final population of the target state Jl), 
does not exceed > 70%. 

In order to design an effective A scheme, i.e. allow¬ 
ing efficient coupling at symmetry, where decoherence 
times are large, we first replace the direct pump pulse 
by a two-photon process, which yields overall the ”2-1-1” 
A scheme [see Fig. 1(a)]. This configuration is however 
known to lack robustness against fluctuations of the pa¬ 
rameters [27, 41-43]. To overcome this problem we sup¬ 
plement the ”2-1-1” A scheme by suitable advanced con¬ 
trol, which turns out to be the key ingredient for achiev¬ 
ing both ~ 100% population transfer efficiency and ro¬ 
bustness. We address two classes of last-generation AAs, 
based on the ” flux-qubit” [2, 44] and on the ’’trans- 
mon” [3, 4, 45] design, respectively. 

We start our analysis from the full Hamiltonian 

H := Ho + Hc{t) (2) 

where Ho := models the undriven AA. The 

control He = QA{t) is operated by a three-tone held 



FIG. 1. (color online) (a) A three-level system with splittings 
uJij := Ej — Ei driven by two quasi resonant ac pump (Sip) 
and Stokes (Sis) fields, in the usual A scheme (ujp := tJo 2 — (5p, 
ujs = Lij \2 — (5s). In 2-1-1 STIRAP the pump is operated by two 
pulses Slpfc at frequencies (Upi := woi — <52, (ap 2 = wi 2 — 5^ + 52 
such that oipi -I- U!p2 = ((^02 — 5 p- (b) Pulses in conventional 
STIRAP (dashed lines) in the counterintuitive sequence, i.e. 
the Stokes pulse is shined before the pump pulse. Real part of 
the pulses in 2-l-l-STIRAP (solid lines): here shows 

the phase modulated control of Eq.(6). 

= J2m=pl,p2,s Am{t) COs[uJmt - ((>m(<)]- ^ IS COUplcd 
to the operator Q, corresponding to the electric dipole 
for natural atoms. In AAs it is, for instance, the charge 
operator in the transmon [3, 4] or the loop current in 
the fiux qudit [2, 44]. Symmetries in the Hamiltonian 
Ho imply that matrix elements Qu = Q 02 = 0. Exter¬ 
nal fields have suitable carrier frequencies (see Fig.1(a)) 
and a slowly varying modulation of the phases (/>m(i), 
for m = s,pl. Rabi angular frequencies are defined as 
■= QoiApi(t), Hp2(t) := Qi2Ap2(t), Hs(t) := 
Qi 2 As(t). For simplicity we take i5p = (5^ = 0 and equal 
peak amplitudes for both the Hpfe(t), where k = 1,2, 
considering Gaussian pulses 

Hs(t) = Hoe“(^) ; Hpfe(t) = e“KV) (3) 

We use the delay [9] r = 0.6 T > 0 which implements the 
counterintuitive sequence. 

Our goal is to reproduce Eq.(l) as an effective Hamil¬ 
tonian yielding STIRAP, by properly shaping the control 
We first consider an AA with a highly anhar- 
monic spectrum, W 12 ^ wio, where each transition can 
be selectively addressed. Therefore we can safely neglect 
the strongly off-resonant Api {Ap 2 and As) in (Iji7j2) 
((OjiJjl)), and also perform the RWA. The three-level 
Hamiltonian in the interaction picture reads 

Ho =i{Hpi(t)e-*^-‘l0)(ll + [Hp2(t)e*[^^‘-^-^(‘)l + 

2 (4) 

+ H,(f)e-*^=(‘)]ll)(2l}+h.c. 

If the two pump pulses are strongly dispersive, j^ 2 |/f^r ^ 
1, they implement an effective two-photon jO) j2) 
pulse which does not populate Jl) [9]. In this regime 
we derive an effective Hamiltonian from the Magnus ex¬ 
pansion of time-evolution operator corresponding to Ho 
[39, 47]. The relevant contributions are found up to 
second order, which captures the coarse-grained dynam¬ 
ics averaged over a convenient time scale At such that 
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Ai|(52| ^ 1 but At fir, At/T, At |(^i(t)| <C 1[39]. Then in 
the same rotating frame of Eq.(l) we obtain 


noise is described by a ’’dissipator” ££> in a Master Equa¬ 
tion of the Lindblad form 


Heff —[{4>p2 — 'Ps) — (82 + 2S'i)]|l)(l|-|- 

+ (^p2 +-S'2 — <S'i)| 2 )( 2 |-|- (-5^ 

+ 2 [(^p|0)(2| + ^^s|l)(2|) + h.c.] 

where ftp{t) = —npinp 2 /{ 2 S 2 ) is the two-photon effective 
pump field and Sk{t) = ( 462 ) are dynamical Stark 

shifts. We see that if we define 6 p := (pp 2 -I- S '2 — and 
:= Pp 2 — Ps — {S 2 + 2Si), Eq.(5) reproduces Eq.(l) 
identifying the effective A system. 

We now look for external control yielding STIRAP. 
It is convenient to take equal pulse amplitudes in H^ff, 
thereby flo = ^r/( 2 |^ 2 |), and the necessary condition for 
adiabaticity [35] sets the time scale T > lO/flo- We fi¬ 
nally adjust the system at both single and two-photon 
effective resonance by choosing the phase modulation ac¬ 
cording to 


PPp2=Si-S2 ; Ps= -(,81 + 282) ( 6 ) 

This is a key point of our analysis: performing the latter 
step is crucial since STIRAP would fail otherwise [see 
Eig.2(a), dashed lines]. Indeed the dynamical Stark shifts 
8 k(t) are of the same order of the effective coupling flp{t). 
Therefore if uncompensated they would determine large 
stray detunings, in particular S{t) = —(82 + 28i) would 
destroy the dark state. 

The phase modulation in Eq.(6) is obtained in closed 
form as a function of the pulse envelopes Am(t) by a 
simple integration. Inserted in the control of the full 
Hamiltonian Eq.(2) it yields the goal we set, namely ~ 
100% efficiency is recovered [see Fig.2(a), solid lines]. 

An important point is that solutions (pmit) of interest 
are slowly varying, consistent with our assumption. This 
is also clear from in Fig. 1(b), where the modulation of 
the Stokes pulse for equal flpiS, i.e. 81 = 82 , is shown. It 
is worth stressing the remarkable agreement between the 
full dynamics and the approximation by H^f f (gray lines 
in Fig.2), which we will use later to estimate appropriate 
figures for Hg and T. 

Noise sources coupled via the operator Q are usu¬ 
ally the most detrimental for decoherence. Effects of 
low-frequency noise from this ’’port” can be suppressed 
by designing a Hamiltonian with suitable symmetries, a 
strategy that has yielded very large decoherence times 
in last-generation superconducting qubits. On the other 
hand high-frequency fluctuations from the Q-port are 
the relevant sources of quantum Markovian noise. Pure 
dephasing is due to residual non-Markovian noise from 
sources coupled to operators orthogonal to Q. The im¬ 
pact of noise is studied using a phenomenological pic¬ 
ture [5, 6, 39], accounting for both Markovian and non- 
Markovian relevant noise sources. Markovian quantum 


pit\S 8 ) = -i[H(Si{t)), pit\Sc.)] + £Dp(t\Si) (7) 


whose solution has to be averaged over a stochastic pro¬ 
cess x(t) describing individual realizations of the non- 
Markovian classical noise. For noise with ~ I//“ low- 
frequency spectrum the leading effect is captured by re¬ 
taining only the contribution of quasistatic stray bias 
x(t) —>■ X of the artificial atom [5, 6], with a suitable 
Gaussian distribution. In this picture stray bias deter¬ 
mine fluctuations of energies AEi and of matrix elements 
AQij, which translate respectively in fluctuations of the 
detunings 5 = A{Ei — Eg) and Sp = A(E 2 — Eq) and of 
the Rabi frequency Hg. Only the former turn out to be 
important [7, 9], thereby Eq.(S7) reduces to the structure 
p = i[p, H{d, (5p)] + CdP, where detunings undergo corre¬ 
lated fluctuations (i5, Sp) induced by x, the full dynamics 
emerging from proper averaging. 

In practical cases a single additional port must be con¬ 
sidered, with associated stray bias x. Then fluctuations 
have a simple linear correlation Sp(x) = aS(x), where 
a is determined by the parametric dependence of the 
spectrum on x (See Ref. 39). In this case experiments 
characterizing the qubit dynamics yield all the needed 
statistical properties of {S, 6 p), since the standard devia¬ 
tion of 5 is 0-5 = V2/T^, where 1/T^ := 1/T* - I/(2Ti) 
is the qubit non-Markovian pure dephasing rate [6] and 
Ti the qubit relaxation time. The multilevel dynamics 
is obtained by averaging over a Gaussian distribution, 
p(S) = ( 2 'Ka'l)~^A^-£‘^ the solution p{t\S,aS) of 
Eq.(S7). We use the Markovian dissipator 


£dP 


^([|l)(l|,p]-2|0)(lHl)(0|) + 
A([|2)(2|,p]-2|I)(2H2)(1|) 


( 8 ) 


accounting for the two allowed transitions in the lowest 
three levels. We assume that Cd does not depend explic¬ 
itly on X, and we retain only spontaneous decay, which is 
the only relevant process at low enough temperature [7]. 
The constant k ~ [Q 2 i/Qio]^ 8 (E 2 — Ei)/ 8 (Ei — Eg) 
depends essentially on the design of the device and, in a 
much weaker way, on the power spectrum 8 (u}), which is 
often ohmic at the relevant frequencies [6, 39]. 

In Fig. 2(a) we present results for the four-junctions 
SQUID of Ref. [2]. They show that A-STIRAP with 
~ 100% efficiency is obtained using T ~ 0.12 ps. We sim¬ 
ulate the dynamics for the lowest six states of the full de¬ 
vice Hamiltonian Hq [39], verifying that leakage from the 
three-level subspace is negligible {J 2 j> 3 Pjj < 2 x lO”"*). 
For this device relaxation (Ti = 12 ps [2]) and the 
associated Markovian dephasing are due to flux noise, 
whereas critical current and charge noise determine non- 
Markovian fluctuations, yielding the overall = 2.5 ps. 
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FIG. 2. (color online) Population histories poo (I) (blue) pii(t) (red) and p 22 {t) (green), (a) For a flux AA biased at the 
symmetry point / = 1/2 (spectrum in the inset), with the phase modulation Eq.6. We used 12^/270 = 200MHz and 82 = —5f2r, 
for the two-photon pump, yielding flo = |I2r/252| = 20MHz. Good adiabaticity, Q.oT = 15, is obtained with T = 0.12ps 
and T = 0.6 T. Results refer to the device Ref. 2 and account for leakage and effects of noise. In the absence of phase 
modulation, population histories an in the absence of decoherence (dashed lines) show no population transfer, (b) Same results 
for a transmon (spectrum in the inset) with phase modulation Eq. (13). Here flo = 3.9MHz from Eq.(ll), T = 0.6ps and 
T = 0.6T. For both designs the approximate effective dynamics (gray thin lines above the exact population histories), obtained 
respectively from Eq.(5) and Eq.(lO), reproduces remarkably well the coarse grained time-evolution. 


We find the remarkable > 97% efficiency, which is essen¬ 
tially limited by Ti only. 

We now turn to AAs based on the transmon de¬ 
sign [3, 39, 45]. Successful implementation of A-STIRAP 
in this class of devices would be very important, since 
they display the largest decoherence times observed so 
far [4, 10], and offer the perspective of fabricating highly 
integrated architectures [13, 14], with a rich arena of ap¬ 
plications. These AAs have a nearly harmonic spectrum 
[inset of Fig.4(b)], quantified by a := 1012 — cuqi and 
P := 0723 — UJ 12 for the four lowest energy levels. Val¬ 
ues of |a| ~ |/3| < Wo 1/10 ensure very large decoherence 
times, at the expenses however of limiting selectivity in 
addressing the desired transitions with strong fields. Har- 
monicity is a severe drawback for operating STIRAP and 
indeed the protocol outlined for flux-based AAs would 
fail in the transmon. In order to find the proper effective 
Hamiltonian we must: (a) include selected off-resonant 
terms of the control, relaxing the quasi-resonant approx¬ 
imation; (b) consider explicitly a fourth level |3) since 
it will determine Stark shifts which must be accounted 
for. We neglect the coupling to the cavity used in the 
transmon as a measuring apparatus and at this stage we 
also assume the RWA, so we consider the Hamiltonian 
H = H 3 + H in the interaction picture, with extra terms 

H ={^^pie-*(^=+“)*[Qi2|l)(2| + Q23e-'^*|2)(3|] 

+ (9) 

+ Q23e-*[^‘-‘^"=(‘)'|2)(3|] -fine.} 

The stray H produces non negligible effects due to the 
fact that anharmonicities |a|, |/3| are small and large Apk 
are needed to yield a sufficient effective dispersive pump 


drive. Since As needs not to be large, the corresponding 
terms can be neglected. A convenient choice of param¬ 
eters turns out to be |i52| ^ |q!|, |/3|. In this regime we 
obtain the following three-level effective Hamiltonian in 
the rotating frame 

-^e// = |9)(2| + |1)(2| + h.c.} 

+ -f (j)p2|l)(l| + (0p2 - <^a)|2)(2| 


where the effective pump coupling is now 


flp — 


HplHp2 O: 

4^2 a-\- 82 


( 11 ) 


and the dynamical Stark shifts of level j due to the cou¬ 
pling to level i under the action of the pk field is given 

by 


sU^) ■■= 


Apk (t) Qij 

2 


1 

^ij ^pk 




■ ^pk 


( 12 ) 


This antisymmetric form for i j accounts also for 
Bloch-Siegert shifts, which are however small in all cases 
treated in this work. Notice that Eq.(lO) includes three 
levels since levels i > 2 only yield Stark shifts. We again 
let dp = (5 = 0, thereby in order to obtain large STIRAP 
efficiency we modulate phases according to 


Pp 2 = T.(^jo-s^ 2 ); - 5^) (13) 

k,j k,j 


If we now use this modulation in the full Hamiltonian 
Eq.(2) we recover ^ 100% efficiency [see Fig. 2(b)]. Again 
the full dynamics is remarkably well approximated by 
the Magnus expansion. Results refer to the transmon of 
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FIG. 3. (color online) Parametric robustness of the protocol, 
(a) Effective duration of the protocol T = IS/flo units) vs 
flr/27r (GHz units) for the transmon of Ref. 4 (Eq.(ll), up¬ 
per curve) and for the fiux-qudit of Ref. 2 (limit a 3> <52, lower 
curve), at fixed |(52|/Hr = 5. (b) Efficiency versus stray detun¬ 
ings S, Sp for the transmon design, showing the robustness of 
the protocol Eq.(13). The solid black inner curve encloses the 
region of efficiency > 95%; results of the approximation Heff, 
Eq.(lO) are also reported (dashed curve) which show again 
the remarkable accuracy of the effective theory. The white 
straight line Sp = 2S represents the correlated quasitatic fluc¬ 
tuations of the stray detunings. The red outer curve encloses 
the > 95% efficiency area for r = 2: it is seen that robustness 
along the line further increases using r > 1. (c) Efficiency 
vs T 2 for the transmon, r = 1 (lower line) and r = 2 (upper 
line). 


Ref.4 and account for both leakage and effects of noise. 
Coherence is again essentially limited by Ti = 70 jjs, 
thereby noise has negligible effects, also allowing for mul¬ 
tiple STIRAP-like cycles. 

Notice that Eq.(ll) implies that the effective peak 
saturates to the value —a/ 2 (S 2 /fir)^, for increasing 17^. 
at constant 62 /^r ^ 1 [see Fig.3(a)]. For this reason the 
duration of the protocol for the transmon [T = 0.6 (is in 
Fig.2(b)] is larger than for the flux-based AA. More gen¬ 
erally, shining larger external fields to shorten the pro¬ 
tocol is useful only to some extent [see Fig.3(a)], but in 
devices with the largest coherence times this is not a lim¬ 
itation. 

Robustness of the protocol is a crucial issue, since 
the success of conventional STIRAP lies in the strik¬ 
ing insensitivity to small variations of control parame¬ 
ters. In the early proposal of ”2-1-1” A—STIRAP, lack 
of efficiency due to the stray dynamical Stark-shift was 
cured by using fields with a small static two-photon de¬ 
tuning [41-43], but unfortunately the resulting proto¬ 
col was not robust [27]. Instead our control scheme 
is tailored to guarantee the same robustness of conven¬ 
tional STIRAP. In Fig.3(b) we show sensitivity against 
fluctuations of the detunings of phase modulated STI¬ 
RAP in the trasmon, which is potentially the most un¬ 
favourable case. For the example shown, frequency fluc¬ 


tuations of the microwave fields < 1 MHz still guaran- 
tee > 95 % efficiency. This important result would be 
hardly attainable for natural atoms driven at optical fre¬ 
quencies [27, 41-43], where the available phase control is 
limited. In addition phase modulated 2-1-1 STIRAP is 
naturally resilient to non-Markovian noise inducing slow 
fluctuations [39] of the energy splittings. This corre¬ 
sponds to fluctuating detunings, correlated as Sp = 2S 
in the transmon of Ref. [4]. The part of this line con¬ 
tained in the high efficiency region of the (5, Sp) plane 
corresponds to T 2 ~ 1 /rs, which sets a figure for the 
resilience of the protocol to non Markovian dephasing. 
Quite interestingly a suitably asymmetric drive with ra¬ 
tio r := Max[Hp(t)]/Max[r2s(t)] > 1 enlarges the stabil¬ 
ity region in a way that low-frequency correlated noise 
affecting the device is dynamically decoupled [see also 
Fig.3(c)]. 

In summary we have shown how to design reliable 
multilevel control in A configuration by 2-1-1 STIRAP. 
The key ingredient is a new control scheme which uses 
pulses with suitable slowly-varying modulated phases, 
Fqs.(6,13). We obtained a unique strategy allowing to 
operate with last generation AAs, where symmetries en¬ 
force selection rules preventing a resonant pump field to 
be coupled directly. It can be easily implemented in such 
devices with available microwave electronics [48], yield¬ 
ing ^ 100% efficiency. It is worth stressing that phase 
control is necessary to guarantee the important property 
of robustness to the same level of conventional STIRAP. 

We finally mention that STIRAP has been very re¬ 
cently observed in the so called Ladder configuration [49, 
50], which is more easily implemented in last-generation 
AAs. It involves a two-photon absorption process, 
whereas A-STIRAP implements a coherent absorption- 
emission cycle. This latter is a fundamental building 
block for advanced control in highly integrated architec¬ 
tures, thereby it would have an impact on applications. 
Phase control exalts in a natural way the advantages of 
last-generation superconducting AAs, where it opens new 
perspectives for advanced quantum control. Our work 
may be extended in these directions using optimal con¬ 
trol theory tools. 
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nostro and J.S. Tsai for useful discussions. 


* gfalci@dmfci.unict.it 

[1] M. Nielsen and I. Ghuang, Quantum Computation and 
Quantum Information (Gambridge Univ. Press, Cam¬ 
bridge, 2010). 

[2] N. Timoney, I. Baumgart, M. Johanning, A. F. Varon, 
M. B. Plenio, A. Retzker, and C. Wunderlich, Nature 
476, 185 (2011). 

[3] J. Q. You and F. Nori, Nature 474, 589 (2011). 









6 


[4] B. P. Lanyon, M. Barbieri, M. P. Almeida, T. Jennewein, 
T. C. Ralph, K. J. Resch, G. J. Pryde, A. OBrien, Jeremy 
L. Gilchrist, and A. G. White, Nat. Phys. 5 (2009). 

[5] J. Pachos, Introduction to Topological Quantum Compu¬ 
tation (Gambridge University Press, 2012). 

[6] H. Kimble, Nature 453, 1023 (2008). 

[7] A. D’Arrigo, R. Lo Franco, G. Benenti, E. Paladino, and 
G. Paid, Ann. Phys. 350, 211 (2014). 

[8] A. Orieux, A. D’Arrigo, G. Ferranti, R. Lo Franco, G. Be¬ 
nenti, E. Paladino, G. Falci, F. Sciarrino, and P. Mataloni, 
Sci. Rep. 5 (2015). 

[9] T. D. Ladd, F. Jelezko, R. Laflamme, Y. Nakamura, 

C. Monroe, and J. L. O’Brien, Nature 464, 08812 (2010). 

[10] M. Devoret and R. J. Schoelkopf, Science 339, 1169 
(2013). 

[11] R. J. Schoelkopf and S. M. Girvin, Nature 451, 664 
(2008). 

[12] P. Macha, G. Oelsner, J.-M. Reiner, M. Marthaler, 

S. Andre, G. Schon, U. Hiibner, H.-G. Meyer, E. Il’ichev, 
and A. V. Ustinov, Nat. Commun. 5, 5146 (2014). 

[13] H. R. Mohebbi, O. W. B. Benningshof, 1. A. J. Taminiau, 
G. X. Miao, and D. G. Cory, Jour. Appl. Phys. 115, 
094502 (2014). 

[14] T. Brecht, W. Pfaff, C. Wang, Y. Chu, L. Frunzio, M. H. 
Devoret, and R. J. Schoelkopf, Npj Quantum Information 
2, 16002 EP (2016). 

[15] T. Niemczyk, F. Deppe, H. Huebl, E. P. Menzel, 
F. Hocke, M. J. Schwarz, J. J. Garcia-Ripoll, D. Zueco, 

T. Hummer, E. Solano, A. Marx, and R. Gross, Nat. 
Phys. 6, 772 (2010). 

[16] Y. Nakamura and T. Yamamoto, IEEE Photonics Jour¬ 
nal 5 (2013), 10.1109/JPHOT.2013.2252005. 

[6] E. Paladino, Y. Galperin, G. Falci, and B. Altshuler, Rev. 
Mod. Phys. 86, 361 (2014). 

[2] J. Bylander, S. Gustavsson, F. Yan, F. Yoshihara, 
K. Harrabi, G. Fitch, D. G. Cory, Y. Nakamura, J.-S. 
Tsai, and W. D. Oliver, Nat. Phys. 7, 565 (2011). 

[4] C. Rigetti, J. M. Gambetta, S. Poletto, B. L. T. Plourde, 
J. M. Chow, A. D. Corcoles, J. A. Smolin, S. T. Merkel, 
J. R. Rozen, G. A. Keefe, M. B. Rothwell, M. B. Ketchen, 
and M. Steffen, Phys. Rev. B 86, 100506 (2012). 

[20] M. Stern, G. Catelani, Y. Kubo, C. Grezes, A. Bienfait, 

D. Vion, D. Esteve, and P. Bertet, Phys. Rev. Lett. 113, 
123601 (2014). 

[21] Y.-x. Liu, J. Q. You, L. F. Wei, C. P. Sun, and F. Nori, 
Phys. Rev. Lett. 95, 087001 (2005). 

[22] M. Mariantoni, M. J. Storcz, F. K. Wilhelm, W. D. 
Oliver, A. Emmert, A. Marx, R. Gross, H. Christ, and 

E. Solano, arXiv:cond-mat/0509737v2. 

[23] J. Siewert, T. Braudes, and G. Falci, Opt. Commun. 
264, 435-440 (2006). 

[24] L. F. Wei, J. R. Johansson, L. X. Gen, S. Ashhab, and 

F. Nori, Phys. Rev. Lett. 100, 113601 (2008). 

[25] J. Siewert, T. Braudes, and G. Falci, Phys. Rev. B 79, 
024504 (2009). 

[26] T. Wilk, S. C. Webster, A. Kuhn, and G. Rempe, Science 
317, 488 (2007). 

[27] K. Bergmann, N. V. Vitanov, and B. W. Shore, Jour. 
Chem. Phys. 142, 170901 (2015). 

[28] P. G. Di Stefano, E. Paladino, A. D’Arrigo, and G. Falci, 
Phys. Rev. B 91, 224506 (2015). 

[29] D. Moller, L. B. Madsen, and K. Molmer, Phys. Rev. 
Lett. 100, 170504 (2008). 

[30] C.-P. Yang, S.-I. Chu, and S. Han, Phys. Rev. Lett. 92, 


117902 (2004). 

[31] Z. Kis and E. Paspalakis, Phys. Rev. B 69, 024510 (2004). 

[32] P. Krai, 1. Thanopulos, and M. Shapiro, Rev. Mod. Phys. 
79, 53 (2007). 

[9] N. Vitanov, M. Fleischhauer, B. Shore, and K. Bergmann, 
Adv. in At. Mol. and Opt. Phys. 46, 55 (2001). 

[34] W. R. Kelly, Z. Dutton, J. Schlafer, B. Mookerji, T. A. 
Ohki, J. S. Kline, and D. P. Pappas, Phys. Rev. Lett. 
104, 163601 (2010). 

[35] K. Bergmann, H. Theuer, and B. Shore, Rev. Mod. Phys. 
70, 1003 (1998). 

[36] Z.-L. Xiang, S. Ashhab, J. You, and F. Nori, Rev. Mod. 
Phys. 85, 623 (2013). 

[37] D. Vion, A. Aassime, A. Cottet, P. Joyez, H. Pothier, 
C. Urbina, D. Esteve, and M. H. Devoret, Science 296, 
886 ( 2002 ). 

[5] G. Falci, A. D’Arrigo, A. Mastellone, and E. Paladino, 
Phys. Rev. Lett. 94, 167002 (2005). 

[39] See supplemental material at [url will be inserted by ed¬ 
itor] for further details on magnus expansion, supercon¬ 
ducting devices and models for dynamics and noise. 

[7] G. Falci, A. La Cognata, M. Berritta, A. D’Arrigo, E. Pal¬ 
adino, and B. Spagnolo, Phys. Rev. B 87, 214515 (2013). 

[41] L. P. Yatsenko, S. Guerin, T. Halfmann, K. Bohmer, 

B. W. Shore, and K. Bergmann, Phys. Rev. A 58, 4683 
(1998). 

[42] S. Guerin, L. P. Yatsenko, T. Halfmann, B. W. Shore, 
and K. Bergmann, Phys. Rev. A 58, 4691 (1998). 

[43] K. Bohmer, T. Halfmann, L. P. Yatsenko, B. W. Shore, 
and K. Bergmann, Phys. Rev. A 64, 023404 (2001). 

[44] C. van der Wal, A. ter Haar, F. Wilhelm, R. Schouten, 

C. Harmans, T. Orlando, S. Lloyd, and J. Mooij, Science 
290, 773 (2000). 

[45] A. Wallraff, D. 1. Schuster, A. Blais, L. Frunzio, R.-S. 
Huang, J. Majer, S. Kumar, S. Girvin, and R. Schoelkopf, 
Nature 421, 162 (2004). 

[3] J. Koch, T. Yu, J. Gambetta, A. Houck, D. Schus¬ 
ter, J. Majer, A. Blais, M. Devoret, S. Girvin, and 
R. Schoelkopf, Phys. Rev. A 76, 042319 (2007). 

[47] J. S. Waugh, L. M. Huber, and U. Haeberlen, Phys. Rev. 
Lett. 20, 180 (1968). 

[48] M. Pechal, L. Huthmacher, C. Eichler, S. Zeytinoglu, 
A. Abdumalikov Jr., S. Berger, A. Wallraff, and S. Filipp, 
Phys. Rev. X 4, 041010 (2014). 

[49] K. Kumar, A. Vepsalainen, S. Danilin, and G. Paraoanu, 
Nat. Comm. 7, 10628 (2016). 

[50] H. K. Xu, C. Song, W. Y. Liu, G. M. Xue, F. F. Su, 
H. Deng, Y. Tian, D. N. Zheng, S. Han, Y. P. Zhong, 
H. Wang, Y.-x. Liu, and S. P. Zhao, Nat. Commun. 7, 
11018 (2016). 



1 


Supplementary Material 

EFFECTIVE HAMILTONIAN BY THE MAGNUS EXPANSION 

We consider a time-dependent Hamiltonian H{t). Our goal is to find an effective Hamiltonian H{t) capturing the 
dynamics on a coarse grained scale, defined by the small but finite time interval At. To this end we write 

n 

U{t + At; t) = Te"* ^ 

fc=i 

where in the latter product we consider n —>■ oo time slices 5tk with X]fe=r = At, we define Hk := H(tk) with tfc 
belonging to the fc-th time slice, and keep ordered in time the exponential operators. By repeated application of the 
Campbell-Baker-Hausdorff relation 


= ex^{A + B + - [A, B] + ...} 


we can write U(t + At;t) = e where i?At(t) is given up to second order by 


HAtit) = 


ft+At 


At 


dt'H{t') - 


2At 


dt' / dt''[H{t'),H{t'')] 


(S2) 


(S3) 


We shall see that in our case this expression can be approximated by a At independent one, i.e. H^t{t) ~ H{t). 

The resulting averaged 17(t -|- At; t) = allows to approximate U{t, to) « Te~*-^*o \ In this way, H{t) is 

identihed as an effective Hamiltonian capturing the dynamics in a coarse grained fashion. 

We carried out our calculations in the interaction picture of Eqs. (4,7) of the main text. In our case each component 
of the pump pulse is far detuned from each transition, allowing us to take At(5(t ^ 1, where d* := \ujij — ujpk\. Since 
^r,\4'pk'{t)\A/T <C 6 fj, we are allowed to choose At such that At/T, AtH^, At|0p2.s(t)| 1. Then the effective 

Hamiltonians Eqs. (5,8) of the main text are obtained from Eq. S3 by bringing out of the integrals all the slowly 
varying terms and subsequently neglecting terms of order {SfjAt)~^ or higher, so that At will not appear in H. The 
physical quantities that this procedure yields are Stark shifts in the diagonal elements and amplitudes for two-photon 
processes in the off-diagonal elements of H. 


SUPERCONDUCTING CIRCUITS 

The superconducting circuits considered in the text are depicted in Fig. SI (a) (flux-based device) and (b) (trans- 
mon). For both we will give the undriven Hamiltonian Hq and the coupling operator Q entering the control Hamil¬ 
tonian He ■= QA{t) [See Eq. (2) in main text] with the associated selection rules and symmetries. 


Flux Qubit 

The flux qubit is made out of a SQUID superconducting loop with four junction. Three of them have equal Josephson 
energies Ej and capacitances C, whereas the fourth one is smaller by a factor a. The Hamiltonian exressed in terms 
of the three independent phases := (</?!, V32, T’a) and the associated charges reads 

3 3 „ 3 3 

H = —Ej cos dpj — olEj cos } + 4 [(1 + 2a) ^ nf — 2 a (S4) 

j j j 

where we neglected the parasitic capacitances [SI]. Here Eq = {20 is the single electron charging energy and 

/ = 4)f,/4>o is the external magnetic flux bias 4)f, in units of the flux quantum $o = h/{2e), and we work at the 
symmetry point / = 1/2. Physically the control is a small magnetic flux ^ac{t) added to 4>h, and A{t) = <I)ac(t)/‘f’o- 
By expanding to first order in A we find the coupling operator Q = 27raEj sin (^^ (/^) = 27rEj///c, where / is 
the loop current operator and Ic the critical current of the big junctions. For / = 1/2 the bare Hamiltonian Hq 
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FIG. SI. Schematics for the superconducting devices considered in the simulations, where (pi are superconducting phases 
across the junctions and the canonically conjugated rii are the numbers of extra Cooper pairs in the islands, (a) The flux-based 
device made out a four Josephson junctions SQUID loop. The three bigger junctions have Josephson energy Ej and parallel 
capacitance C while the smaller one has Josephson energy and capacitance smaller by a factor a. (b) The Cooper pair box 
is a superconducting island separated from the circuit by a single Josephson junction of energy E,j and capacitance C. The 
transmon is a capacitively-shunted Cooper pair box with Ej / Ec ^ 1 strongly coupled to an electromagnetic transmission line 
resonator, which allows quantum non-demolition measurements. 

enjoys a symmetry with respect to the parity operator Vip\ip) = |—<p). As a consequence, eigenfunctions |n) of Hq 
can be chosen with a definite symmetry, i.e. <^) = (—!)"■!/'„(<^), where n labels eigenenergies in increasing order, 

implying the selection rule for the odd parity coupling operator Q 

(n|Q|m) oc I-(-1)’"+" (S5) 

In our simulations, we employed values EcI‘2.'k = 4GHz, Ej = bQEc and a = 0.54 as in Ref S2, and used 13^ charge 
states to diagonalize Hq. 


Transmon 

The transmon can be modelled by a Cooper pair box [see Fig. SI (b)] in the Ej/Eq ^ 1 regime [S3]. Here Ej is the 
Josephson energy of the junction of the box, while the charging energy Eq = {20tot involves the total capacitance 

Ctot = C + Cg. The Hamiltonian reads 

H = —Ej cosip + AEcin — qg — A{t))^ (S6) 

where qg = CtotVgl2e is a bias parameter. The control field A{t) = CtotVac{t)I2e couples to operator Q = —8Ec{n — 
qg). Biasing the device at qg = 1/2 the Hamiltonian Hq is symmetric with respect to the charge-parity operator 
Vn = |1 ~ w)(m|, implying that also in this case a selection rule of the kind Eq.(S5). In the simulations we 

employed values EcI2ti = 0.212GHz, Ej/Ec = 49 as in Ref. [S4], using 100 charge states to diagonalize Hq. 


EFFECTS OF NOISE 

Different noise sources act on the devices, but independently on their nature they produces essentially two distinct 
classes of effects [S5, S6]. Environmental modes with frequencies comparable to Bohr or Rabi frequencies act as 
sources of Markovian quantum noise, whose leading low-temperature effects in STIRAP are spontaneous decay, the 
associated secular dephasing and field-induced absorption [S7, S8]. At lower frequencies noise in the solid state is 
non-Markovian and exhibits a l//“ behavior. The leading effect is a pure dephasing, analogous to inhomogeneous 
broadening produced by a classical noise source, and it is effectively described by a stochastic external drive with the 
desired spectral features [S5, S6]. 

According to this picture, we describe Markovian quantum noise by a standard dissipator Cd term in a Lindblad 
Master Equation 

P/(i|x) =-i[H{SiL{t)),pf{t\SiL)\+CDPf{t\SiL) (S7) 

where we also account for the effect of non-Markovian noise by allowing the parametric dependence on the classical 
stochastic processes {x(t)}. Averaging over this latter yields the full noisy dynamics 

P{i) = J P[i(t)] p{t\Si) 


(S8) 
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FIG. S2. Bohr frequencies Ei — Eq (lower coloured curve) and E 2 — Eq vs Ej/Ec (upper coloured curve) around a bias point J 
for a flux qudit (a) and a transmon (b). It is seen that Bohr frequencies vary linearly with respect to the fluctuating parameter 
Ej/Ec- Variances us and as^ are indicated in correspondence of a noise variance aj. Correlation of detuning fluctuations, i.e. 
Sp = aS, is determined through a — [d{E 2 — Eq) /d{Ej/Ec)]/[d{Ei — Eq) /d{Ej/Ec)], where the derivative is evaluated at the 
bias point J. 


If low-frequency noise has a ~ l//“ spectrum the main effect is to induce quasistatic stray bias x of the artificial 
atom. The path integral Eq.(S8) can therefore be evaluated in the Static Path Approximation(SPA) [S5, S6], reducing 
to an ordinary integration over random variables x. These latter are moreover Gaussian distributed if we assume that 
they are due to many uncorrelated microscopic sources. 

We apply this recipe to 2-1-1 STIRAP in HNP devices. The emerging important qualitative issue is that noisy 
three-level dynamics is fully characterized by decoherence in the ’’trapped” (or qubit) subspace only, plus information 
on the Hamiltonian of the device alone, a results also obtained conventional STIRAP [S7]. Indeed stray bias due to 
low-frequency noise determine fluctuations of energies and of matrix elements AQ^-. They translate respectively 
in fluctuations of the detunings 5 = AEi — AEg and 5p = IS.E 2 — AEg and in fluctuations of the Rabi frequency Hg. 
The sensitivity to such parameters has been extensively studied [S9]: while fluctuations Hg are irrelevant for STIRAP, 
fluctuations of detunings are important. Therefore the relevant open system dynamics turns out to be described by a 
Lindblad Master Equation with the structure p = i[p,H(S, Sp)]+CdP where H depends on fluctuations {S, Sp) induced 
by stray bias x. In Fig.3(c) of the main text we plot the efficiency of the protocol vs stray (<5, Sp): efficiency is large 
if fluctuations do not let the system diffuse out of the central diamond region. In particular the protocol is critically 
sensitive to fluctuations of the two-photon detuning 6, i.e. of the ’’qubit” splitting Ei — Eq. 

Concerning quantum noise, CdP includes in principle the various decay rates and associated excitation and secular 
dephasing processes, but in practice again the ’’qubit” spontaneous decay only has to be accounted for, i.e. the only 
relevant term turns out to be CdP = p] ~ 2cr_pcr_|_) where 1/Ti is the spontaneous decay rate of the 

qubit, and the Lindblad operators are the corresponding lowering and raising operators tT_ = = |0)(1|. Indeed 

selection rules suppress 2^0 processes, whereas we can estimate the transitions rate r 2 _>.i = k/Ti, where k = 0(1) 
depends on features of the device, and directly check that in the devices of interest it has no impact on the results. 
This is due to the fact that |2) is depopulated when STIRAP is successful. The same holds true for field induced 
absorption, since {^pk) act when |1) (|0)) are depopulated. Notice finally that we did not include in the dissipator 
extra pure dephasing terms since they are accounted for by the average over non-Markovian quasistatic fluctuations. 

We now briefly discuss how figures of noise can be extracted from experiments on qubits, referring to the flux qudit 
of Ref. S2. The dominant source of decoherence is flux noise. The low-frequency part of its spectrum 5 '$(cli), though, 
has minimal effect since at the symmetry point / = 1/2, energy fluctuations are quadratic in the small corresponding 
stray bias xi/t) = f(t). Instead its high-frequency components determine the qubit Ti = 12 ps and the rate r2_>.i. For 
this latter we estimated k = [Q 2 i/Qio]^ [«S'$[(E 2 — Ei)/Ii]/S'$[(£’i — Eo/jK) < 1, taking for S'$(w) the linear behavior 
observed for quantum noise [S2], and verified that it does not affect the dynamics. Subdominant noise sources for 
flux qudits are critical current and charge noise. They do not produce relaxation at the symmetry point, but they are 
the main source of non-Markovian pure dephasing. The induced energy fluctuations are linear in the corresponding 
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stray bias x(t) (see Fig.S2), determining Gaussian suppression of the qubit coherence [S 6 ]. The power spectrum 
5'x(w) ~ 1/w has been measured from the qubit dynamics, over several frequency decades [S2]. Therefore we can 
safely use the SPA result for the qubit coherences decay oc and extract from the measured non-Markovian 

qubit pure dephasing rate I/T 2 := l/TJ ~ l/(2Ti) the variance as = of the two-photon detuning in STIRAP. 

Fluctuations of 6 p, i.e. of E 2 — Eq, are easily found from the parametric dependence on the external bias of the 
calculated spectrum of the device (see Fig.S2). Notice that since each source of noise induces a single stray bias Xi, 
fluctuations of detunings are correlated. For HNP devices subdominant noise sources induce fluctuations Sp = aS. 
The constant a depends on the band structure of the device (see Fig. S2) and in particular a = —4.5 in the flux qudit 
of Ref. S2 refers to critical current and charge noise. 

Similar considerations hold for the transmon of Ref. [S4] a ~ 2 and fc ~ 2, symmetry suppresses low-frequency 
charge noise, and subdominant noise as flux and critical current noise lead to pure dephasing I/T 2 which is very 
small. 

Numerics have been carried out through a Montecarlo quantum jump approach accounting for Markovian noise, by 
averaging over 10^ trajectories. Non-Markovian noise has been taken into account by a further average: we impose 
Sp = aS and sample S from its Gaussian distribution for each trajectory. Variances we used are (T 5 = 4.1 x 10“^ flo = 
—0.22 as for flux qubit (flo = 20 MHz) and = 2 x 10“"^ Oo = 0.5 cr^ in the transmon. By inspection of Fig.3(c) 
of the main text it is clear that STIRAP in HNP devices is robust against such low.frequency fluctuations for both 
devices, as it also results from simulations. 
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